Spin relaxation and decoherence of holes in quantum dots 



in 
o 
o 



oo 



c3 



I 

C 

O 

o 



> 

00 

m 
O 
in 
o 

-i— > 
a 



i 

c 

o 
o 



Denis V. Bulaev and Daniel Loss 
Department of Physics and Astronomy, University of Basel, 
Klmgelbergstrasse 82, CH-4056 Basel, Switzerland 
(Dated: February 2, 2008) 

We investigate heavy-hole spin relaxation and decoherence in quantum dots in perpendicular 
magnetic fields. We show that at low temperatures the spin decoherence time is two times longer 
than the spin relaxation time. We find that the spin relaxation time for heavy holes can be com- 
parable to or even longer than that for electrons in strongly two-dimensional quantum dots. We 
discuss the difference in the magnetic-field dependence of the spin relaxation rate due to Rashba 
or Dresselhaus spin-orbit coupling for systems with positive (i.e., GaAs quantum dots) or negative 
(i.e., InAs quantum dots) p-factor. 



Spin physics has become one of the most rapidly de- 
veloping branches of condensed matter physics. Spin 
physics is very important, not only from a fundamental 
point of view, but also for the fabrication of novel elec- 
tronic devices, for the experimental realization of quan- 
tum computation, and for the development of spin elec- 
tronics (spintronics) 0. Quantum dots (QDs) are most 
attractive candidates for these applications because of 
their reduced dimensionality, leading to long-lived spin 
states and allowing single spin manipulation |2( . 

Recent experiments 0, El 01 show that electrons in 
QDs have a long spin relaxation time (up to 20 ms 
and it is now possible to prepare a single electron spin 
state with a well-defined orientation, read the spin state 
out, and store the information about the spin orientation 
for a long time 0]. There are two main spin relaxation 
mechanisms for electron spins in QDs: that due to the 
electron-phonon interaction @, IS E| an d that due to 
the hyperfine interaction with surrounding nuclear spins 
U3 EH E3- Since the valence band has p symmetry, 
the hyperfine interaction of holes with lattice nuclei is 
suppressed with respect to that of the conduction band 
(electrons). This has led to an increased interest in hole 
spins as carries of long-lived quantum information. It was 
shown that in thin quantum wells (QWs) the hole spin re- 
laxation is slower than that in the bulk case [T^Hil ]. Nev- 
ertheless, the hole spin relaxation time is several orders 
of magnitude smaller than that for electrons. This is due 
to the fact that, in addition to existing spin-orbit (SO) 
couplings for electrons due to bulk inversion as ym metry 
(BIA) (Dresselhaus spin-orbit (DSO) coupling |lj|) and 
structure inversion asymmetry (SIA) (the Rashba spin- 
orbit (RSO) coupling pjj) there is strong SO coupling 
between the heavy-hole (HH) and light-hole (LH) sub- 
bands [13 • 

Very recently, investi gatio n of hole spin relaxation in 
QDs was reported [Isl llflj. In these works only one 
SO mechanism was considered, the SO coupling between 
HHs and LHs. It was shown that the hole spin relaxation 
time in QDs is longer than that in QWs but still shorter 
by several orders of magnitude than that for electrons 
in QDs. Furthermore, it was found that SO coupling 



between HHs and LHs is negligible for two-dimensional 
(2D) QDs if the energy splitting between the HH and LH 
subbands is much larger than the level spacing in those 
subbands [1^, Up to now this case has not been investi- 
gated, though it is very important for the realization of 
coherent spin states with a spin relaxation time longer 
than that for electrons. In this case other SO coupling 
mechanisms (RSO and DSO couplings) become signifi- 
cant. 

In this Letter, we study HHs confined to a QD in a per- 
pendicular magnetic field. We consider the three main 
SO coupling mechanisms: RSO, DSO, and SO coupling 
between the HH and LH subbands. An effective Hamilto- 
nian for 2D HHs is derived. We study the spin relaxation 
and decoherence of the system induced by the interaction 
of HHs with phonons. 

From the two-band Kane model, the Hamiltonian for 
the valence band of III-V semiconductors is given by 
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where Hlk is the Luttinger - Kohn Hamiltonian |2lJ, 7 
is due to BIA, 77 = A so /(E g + A so ), A so is the split-off 
gap energy, E g is the band gap energy, J = (J x ,J y , J z ) 
are 4x4 matrices corresponding to spin 3/2, Q, z — 
P z {Px — Py )> and ^x, &y are given by cyclic permuta- 
tions. The last term in Eq. QJ, caused by SO interaction 
of the conduction and valence bands, is DSO coupling 
for the valence band (for the conduction band, the DSO 
coupling is given by |lj| 7<r • ft, where cr — (a x , a y ,a z ) 
is the vector of Pauli matrices). The magnetic field 
induces a Zeeman splitting, which is described by the 
following term 22]: H z = -2K/i fl B • J - 2 9 /xrB • J, 
where k and q are the Luttinger parameters [22J and 
3 = ( J x> J y> J z)' For 2D asymmetric QWs, due to 
SIA along the growth direction, there is an additional 
SO term, the RSO term, which in the two-band model 
is given by [H 13 a R P x E • J, where a R is the RSO 
coupling constant and E is an effective electric field along 
the growth direction. 

We consider a [001]-grown 2D system. Due to con- 
finement along the growth direction, the valence band 
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splits into a HH subband with J z — ±3/2 and a LH sub- 
band with J z = ±1/2 17]. If the splitting of HH and 
LH subbands is large, we describe the properties of HHs 
and LHs separately, using only the 2x2 submatrices for 
the J z — ±3/2 and J z — ±1/2 states, respectively. The 
HH submatrices have the property that J x = J y = 



and J z — |cr z (see Ref. I25TI . For such a system and 
low temperatures only the lowest HH subband is signifi- 
cantly occupied. In this case, we consider HHs only. In 
the framework of perturbation theory , using Eq. flj 
and taking into account the Zeeman energy, the RSO, 
and DSO term, the effective Hamiltonian for HHs of a 
QD with lateral confinement potential U(x,y) is given 
by 
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H = ^APt + P$) + U(x, y) + H% - -gfi B B ■ <x, (2) 

where m — too/(7o + 7i) is the effective HH mass, tuq 
is the free electron mass, 70 and 71 are the Lutinger pa- 
rameters [z| , g is the effective g-factor of HHs, and 
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is the SO coupling of HHs, which is due to the SO cou- 
pling between LH and HH subbands, SIA (the first term), 
and BIA (the second term). The first term in Eq. © 
is the RSO coupling jU 0] and the second term is 
the DSO coupling. Here, a = 3^oa.R(E z ) /2m A, (3 = 
37o7(-P?)/2too^A, <7± = (a x ± ia y )/2, P ± = P x ± iP y , 
(E z ) is the averaged effective electric field along the 
growth direction of a QD, and A is the splitting of LH 
and HH subbands. The splitting between HH and LH 
subbands A ~ d~ 2 , where d is the QD height. Compar- 
ing Eq. with the SO coupling term for electrons, we 
find that for a QD with the characteristic lateral size I, 
the ratio (H? )/{H%) ~ (l/d) 2 . Thus, for strongly 2D 
QDs (I ^> d), the SO coupling of HHs can be less than 
that for electrons. 

Without the SO interaction (a = (3 = 0), the spectrum 
of a QD Hamiltonian with parabolic lateral confinement 
can be found using a canonical transformation of the 
Hamiltonian |26| . For circularly-symmetric QDs with 
oscillator frequency uj (U(x,y) — raus^ix 1 + y 2 )) in a 
perpendicular magnetic field, the energy spectrum and 



huj 



wave functions of HHs are given by |8j E, 
hw-{nx + 1/2) 

m, ri2 



0,1,2, 



(n 2 + 1/2)t fiwz/2, \nm 2 T (I)) = 
n2 (q 2 ^mu + /h) |t (I)), where 
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£1 ± uj c /2 are the hybrid 
= \e\B/mc is the 



frequencies, ft = y/tjQ + uf/i, 
cyclotron frequency, loz = g^BB/h is the Zeeman fre- 
quency, q\ and go are new coordinates in the transformed 
phase space |26|, and <& n (q) are oscillator functions. 
We only consider low-lying levels. For definitcness, wc 
assume that B > and g > 0, then the ground state 
is the spin-up state 1 00 f). At low B, the next level is 
Eqoi, which is Zeeman-split from the ground state level, 



and at high B, levels E n0 ^ are close to the ground state 
level (since w_ — > as B — > 00). Therefore, there are 
crossings of levels E n0 ^ with E Q0 ^ at u>z = nu>_. 

We now take the SO coupling of HHs into account. 
The spectral problem for H can be solved in the frame- 
work of perturbation theory The SO coupling influ- 
ences the wave functions more than the energy spectrum 
(since the energy corrections due to are only second 
order). Thus, SO coupling leads to mixing of spin-up 
and spin-down states. The RSO and DSO terms dif- 
fer by symmetry in momentum space 0, 127I ] and hence 
mix different states: the state |00 |) mixes with |03 J.) 
(E031 — -E-oot = ft(3cj + + u>z)) due to RSO and with the 
states |01 |) {E m - E m = h{uj + + u z )) and |12 |) 
due to DSO coupling. In turn, |00 |) mixes with 1 30 f) 
(E 301 -E oi = h{3uj^-ujz)) due to RSO and with 1 10 j) 
(E101 — E oi = H(lu- — uz)) and |21 f) due to DSO cou- 
pling. Again, we only consider the case B > 0, since the 
physical properties of the system are independent of the 
sign of B 8J. In this case, u + > \u>z\ and the mixed state 
levels cross (at loz = (J- and luz = 3w_) for positive HH 
g-factor (ojz > 0) but not for g < 0. Therefore, for 
g > (e.g., GaAs QDs |25|1. there is strong spin mixing 
of the states at these points and the SO term J3J leads to 
anticrossings of the corresponding levels ||| ■ Strong mix- 
ing of spin-up and spin-down states and an anticrossing 
at uz = (J- ((Jz = 3cj_) are due to only DSO (RSO) 
coupling. For g < (e.g., InAs QDs |28|), the levels 
i?oot and Eiq± (-Eoot and U304.) just cross each other at 
uz = (J- ((Jz = 3 J- ) , since there is no coupling between 
the corresponding states. 

In the following we study spin relaxation induced by 
phonon-HH interactions only. The coupling between HHs 
and phonons with mode qa (q is the phonon wave vector, 
and the branch index a = L, Tl, T2 for one longitudinal 
and two transverse modes) is given by l2flj | 



K h a (r) 
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where p is the crystal mass density, s a is the sound 
velocity, V is the volume of the QD, q|| = (q x ,q y ), a 
and b are the constants of the deformation potential, 
Aqa = ^i£,jdf a f3iji, £ = q/<7, d qQ is the phonon polariza- 
tion vector, and fliji is the piezotensor, which has nonzero 
components only when all three indices i, j, I are different: 
fixyz = flxzy — ■ ■ ■ — hu/es (ss is the static dielectric 
constant). For GaAs, ehu = 1.2 x 10 7 eV/cm, £5 = 13.2, 
and for InAs, ehu = 3.38 x 10 6 eV/cm, and e s = 14.6. 
In Eq. we introduced the form- factor F(q z ) which is 
determined by the spread of the HH wave function in 
the z-direction: F(q z ) = J dze iqzZ \ipo( z )\ 2 , where tpo(z) 
is the ground state envelope wave function of a HH along 
the z-direction. The form factor F(q z ) equals unity for 
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\q z \ <C d^ 1 and vanishes for \q z \ » d~ x 0. 

We consider a single-particle QD, in which a HH can 
occupy one of the low-lying levels. As mentioned above, 
with increasing B some energy levels with the same spin 
orientation cross the upper Zeeman-split ground state 
level and we should study the relaxation of an n-level 
system, the first n — 1 levels have the same spin and the 
n-th level has the opposite spin orientation. In the frame- 
work of Bloch - Redfield theory ji^] , the Bloch equations 
of HH spin motion for such a system in the interaction 
representation are given by 

(S z ) = (Sr-iS^/Tx-Rit), (5) 
(S x ) = -(S x )/T 2 , (S y ) = -(S y )/T 2 , (6) 

where R(t) = W n ip nn {t) + J27=i w „ipu(t), p(t) is the 
density matrix, Wij is the transition rate from state j to 
state i, St is a constant (which hase the value of (S z ) in 
thermodynamic equilibrium if R(t) = 0), 
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1 
2J\ 
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^£^i, (7) 
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where the pure dephasing (due to fluctuations along z 
direction) is absent in the spin decoherence time T 2 since 
the spectral function is superohmic. As can be seen from 
Eq. JS}, the spin motion has a complex dependence on 
the density matrix and, in the general case, there are 
n — 1 spin relaxation rates. However, in the case of low 
temperatures (hqs a 3> T), when the phonon absorption 
becomes strongly suppressed, solving the master equa- 
tion, we find that R(t) w 0, therefore, there is only one 
spin relaxation time T\. In this limit, the last sum in 
Eq. is negligible and the spin decoherence time satu- 
rates, i.e., T 2 = 27\. 

For brevity, we present only the probability W\ n of 
transition with phonon emission for the Zeeman-split 
ground state which can be expressed as a sum of two 
terms due to RSO and DSO couplings, respectively: 
Wm = + W? w where 

2 



W-r 

In 



wP 

vv In 



2 h\ 



AN ug + 1) 



x E s « 9e ^ |i2/2s ' /(7) ' 



3cj- — ujz 



f] 2 m 2 huj 

f2j(3) 



1) 



£• 



2s, 



a + ~ \ £ ■ d qQ — ht z d* a 



In the case of parabolic confinement along the growth 
direction of a QD, can be expressed in terms of 
error functions 8]. As mentioned above, for g > 0, 
the SO term leads to level anticrossings at u>z = 
iv- (due to DSO) and ojz — 3w_ (due to RSO). In 
this case, the denominators 3w_ — ujz in Eq. (J8J and 
cij_ — u>z in the expression for / should be replaced by 
sgn(3w_ — luz)\/{3u- — wz) 2 + (Ar/H) 2 and sgn(w_ — 
loz)\J{u- — ujz) 2 + (A_d/?i) 2 , respectively. Here = 
2/3(mZ) 3 w_(aj 2 + w 2 ) and = 2\/6a(mluj^) 3 are the 
level splittings at the anticrossing points. 

Note that there is no interplay between RSO and DSO 
couplings for HHs of a QD in perpendicular magnetic 
fields, as is true also for electrons 

Figure ^ shows the dependence of the spin relaxation 
rate for HHs in a GaAs (see Fig.^( a )) an d InAs (see 
Fig. 12(b)) QD on a perpendicular magnetic field. The 
solid (dashed) curve corresponds to the spin relaxation 
due to RSO (DSO) coupling. In the case of a positive HH 
(/-factor (i.e., a GaAs QD), there are peaks (see Fig.U(a)) 
in the relaxation rate curve at loz — w_ (ujz — 3w_), 
which are caused by strong spin mixing at the anticross- 
ing points due to DSO (RSO) coupling Q. The half- 
width of the first peak is Ad/H and that of the second 
peak is Ar/H. Since, at anticrossing points, ui + ^> 
the half-width of the peak, which is caused by the DSO 
coupling, is much larger than that due to the RSO cou- 
pling (in our case, A# = 4 neV and Ad = 1.5 /ieV). In 
the case of a negative HH g-factor (i.e., an InAs QD), 
there is no mixing of spin-up and spin-down states of 
the crossing levels, therefore, the field dependence of the 
relaxation rate is monotonic (see Fig.^(b)). 

From Eqs. JHJ and it can be shown that at low 
temperatures and at low magnetic fields (B < 0.5 T) 
the relaxation rate due to RSO (DSO) coupling ~ B 9 
(~ B 5 ). Therefore, the field dependence of T\ due to 
DSO coupling of HHs in a QD is the same as for electrons 
0001 but that due to RSO coupling is quite different. 
This qualitative difference in the field dependence of the 
(8) spin relaxation can serve to provide information about 
the leading SO interaction term at low B. 

Let us consider the quantitative difference between the 
spin relaxation of electrons and HHs. For simplicity, we 
consider DSO coupling only. It can be shown that at low 
B (B < 0.1 T) the ratio of the HH spin relaxation time 



.(9) 



and that for electrons T-f is given by 



Here N u = (e hul l T - l) , I = y/h/mSl, f = {uj 2 _ + 
w\)[uj + I '(u++u>z)-u- 1 '(u>--ljz)], j = lo-lo + [uj + /(lo- + 
2uj+ + ujz) - w_/ (2uj- + uj + - ujz)], 
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where g e i and m e i are the electron (/-factor and effective 
mass, respectively. Usually the (/-factor and the effective 
mass of an electron are less than those of a HH and the 



4 



10' 



10' 



10" 



10" 



~ 10" 

t 

~ 10 3 



10 
10 

1 




r (b) 

10 
















r ^ 1 






r 0.1 










= 0.01 0.1 1 10 




B(T) 






















7 I 




I I 



10 



B(T) 



FIG. 1: (color online). Spin relaxation rate 1/Ti of HHs (solid and dashed lines) and electrons (dotted line) in a GaAs (a) and 
an InAs (b) QD 



= 5 nm, lo = \J li/rnujo — 30 nm, and T 
(a = P). For a GaAs QD (a), g = 2.5 [H, 7/ft 3 = 28 eV A 3 
7/ft 3 = 130 eV A 3 Uni, and A = 150 meV. We note that T 2 = 
relaxation time due to DSO coupling. 



= 0.1 K) due to DSO (solid line) and RSO (dashed line) coupling 
23, and A = 40 meV [HEl, for an InAs QD (b), g = -2.2 0, 
= 2Ti. Insets: Ratio between the HH (Ti) and electron (Tf) spin 



spin relaxation time Tf l for electrons is much longer than 
for HHs. However, for strongly 2D QDs (Iq ^> d), the 
spin relaxation time for HHs can be comparable to, or 
even longer than, that for electrons. For the GaAs QD 
considered here (77 = 0.18), T\ for HHs is comparable 
to that for electrons 0, 0] (see inset Fig. H] (a)). Now 
consider the InAs QD, for which g e i is larger than for the 
GaAs QD and the energy gap is narrow (77 = 0.48). We 
find that the spin relaxation time for HHs is longer than 
for electrons (see inset in Fig.^(b)). 

Since a and (3 ~ d 2 , the spin relaxation rate increases 
with increasing confinement along the growth direction 
of a QD and, as can be seen from Eqs. JSJ and ©, with 
decreasing lateral confinement (with a decrease in the 
confinement frequency loq). As follows from Eq. J2J, the 
HH spin relaxes primarily due to piezoelectric phonons 
at low B and due to deformational acoustic phonons at 
B > IT. 

In conclusion, we have shown that due to the different 
symmetries of the RSO and DSO terms in momentum 
space, these terms lead to different behavior of the spin 
relaxation: at low magnetic fields T\ ~ B~ 9 in the case of 
RSO coupling and T\ ~ B~ 5 in the case of DSO coupling. 
The field dependence of the spin relaxation rate is mono- 
tonic for a system with a negative g- factor (i.e., HHs in 
InAs QDs), whereas for g > (i.e., HHs in GaAs QDs) 
the relaxation rate has peaks corresponding to level anti- 
crossings and the associated enhanced mixing of spin-up 
and spin-down states at the anticrossings. 
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